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ABSTRACT 


Under  the  assumption  of  constant  failure  rates  it  is  pos¬ 
sible  to  build  a  "reliability  shorthand"  which  gives  a  sim¬ 
ple,  unified  approach  to  reliability  computations  for  sys¬ 
tems  in  the  presence  of  complications  like  support  by  shared 
spares  or  changes  in  the  failure  rates  of  surviving  compo¬ 
nents  when  other  components  fail.  The  computational  imple¬ 
mentation  of  the  shorthand  depends  upon  the  convolution  of 
strings  of  exponentially  distributed  random  variables. 

This  paper  presents  an  algorithm  for  the  numerical  convo¬ 
lution  of  exponentially  distributed  random  variables.  After 
reducing  the  system  scenario  to  its  shorthand  format,  one 
can  use  the  programs  that  are  given  in  the  appendix  to 
obtain  numerical  values  for  the  reliability  of  the  system. 
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I.  IMSQB2 CIIQS 


The  reliability  shorthand  considered  in  this  paper  has 
been  developed  in  conjunction  with  the  course  OA  4302  "Reli 
ability  And  Weapon  Systems  Effectiveness  Measurement"  at  th 
Naval  Postgraduate  School.  A  tutorial  introduction  to  the 
reliability  shorthand  was  given  by  Repicky[  Ref .  2  ].  This 
paper  is  devoted  to  a  complementary  part  of  the  idea. 

Any  study  on  system  reliability  always  requires  two 
steps;  one  is  the  description  of  the  system’s  life  and  the 
other  is  the  derivation  of  its  survival  function. 

Under  the  assumption  of  constant  component  failure  rates 
this  paper  presents  a  way  of  obtaining  the  system's  reli¬ 
ability  which  requires  little  beyond  the  description  of  the 
system's  life. 

Section  II  contains  an  approach  to  the  convolution  of 
exponentially  distributed  random  variables.  Also  there  is 
presentation  of  a  computational  algorithm  for  the  convolu¬ 
tion  of  exponentially  distributed  random  variables. 


Appendix  A  gives  survival  functions  corresponding  to  sev¬ 
eral  reliability  shorthand  notations  and  a  program  writen  in 
Fortran  for  computations  from  shorthand  notations. 

Section  III  deals  with  the  reliability  of  redundant  sys¬ 
tems  under  the  assumption  of  constant  component  failure 
rates. 

Appendix  B  consists  of  a  program  written  in  Fortran.  The 
program  supports  the  approach  of  Section  III.  There  is  a 
crude  Bonte  Carlo  simulation  program  in  Appendix  D  which  is 
a  simulation  program  parallel  to  the  program  in  Appendix  B. 

There  is  another  program  in  Appendix  C  written  in  For¬ 
tran.  This  program  uses  the  network  approach  to  systems 
described  in  Section  III. 

Appendix  B  summarizes  the  definitions  used  in  this  paper. 


II. 


||^ A^|HO AC|^|C^CO MP g £ING  CONVOLUTIONS  OF  EXPONENTIAL 


This  section  introduces  a  general  algorithm  for  com¬ 
puting  the  survival  function  of  any  convolution  of  exponen¬ 
tial  random  variables. 

In  reliability,  the  teem  convolution  refers  to  the  summa¬ 
tion  of  independent  random  lives.  In  order  to  have  simplic¬ 
ity  in  specifying  convolutions,  the  reliability  shorthand 
introduces  a  special  notation.  In  the  following  sections  we 
will  use  this  notation. 

A.  THE  SURVIVAL  FUNCTION  POR  A  CONVOLUTION  OF  RANDOM 

VARIABLES 

• 

Let  FI  (t)  and  F2(t)  be  the  survival  functions  for  the 
random  variables  T1  and  T2  respectively.  Let  fl(t)  and  f 2  (t) 
be  the  corresponding  densities.  Let  P  (t)  be  the  survival 
function  for  the  random  variable  T,  where  T=T1+T2. 

Then  the  likelihood  expression  for  F  (t)  is 
t 

F(t)=F1  (t)  *f?2  (t-s)  ?1(s)  ds 


E 


In  the  right  hand  side  of  equation,  FI  (t)  is  the  prob¬ 
ability  that  component  one  completes  the  mission, 
t 

Ji~? 2 (t-s)  f  1  (s)  ds  is  the  probability  that  at  some  time 
s  (0£s<t)  coaponent  one  fails,  coaponent  two  takes  its  place 
and  carries  the  system  to  the  end  of  the  mission  duration  t. 

In  order  to  illustrate  consider  some  applications. 

1 .  Example 

Beliability  shorthand  Notation  :  EXP  }  +EXP  {X2 } 
SYSTEM:  One  coaponent  having  one  spare  with  a  dissimilar 
failure  rate.  If  the  active  component  fails,  the 
spare  will  replace  it  immediately. 

Here  the  life  for  the  3ystea  is  T  =  T1  +  T2.  The  reli¬ 
ability  shorthand  notation  indicates  that  this  system  has  an 
exponential  life  with  failure  rate  Aa  followed  by  an  expo¬ 
nential  life  with  failure  rate  \z . 

The  survival  function  for  the  active  component  is 

F-,  (t)=  e'Xlt  ,  t>0. 

The  survival  function  for  the  spare  is 
\  (t)  =  e‘Xzt  ,  t>0. 

The  survival  function  for  the  system  is 


The  result  comes  out  as  expected  to  be  the  Erlang  (2,X  } 
survival  function. 


B.  RELIABILITY  SHORTHAND  NOTATION 

EXP(^}  +  EXP  {!}♦.... +EXP{X}  (n  identical 

exponential  lives) 

This  leads  to  the  Erlang  {n,  X  }  survival  function 
n 

F(t)  =  l  (Xt)  /(i-1)  !  e'  ,  t>D. 

1=1 

The  use  of  moment  generating  function  gives  the  result 
immediately. 

C.  THE  RELIABILITY  SHORTHAND  NOTATION 

EXP  {Ai }  +EXP  (  X2}  ♦ . +  EXP  (Xn) 

This  is  the  expression  for  the  convolution  of  n  random 
variables  where  each  has  a  distinct  failure  rate. 

By  the  approach  used  in  Section  2.1,  adding  one  exponen¬ 
tial  life  at  once,  one  can  obtain  zhs  formula  for  the  survi 


val  function 


This  is  also  a  well  known  formula  and  can  be  obtained 


from  by  a  moment  generating  function. 

D.  THE  RELIABILITY  SHORTHAND  NOTATION 

EXP{Xi)  ♦EXP{X2}  + . ♦EXP{X2}  (n4  terms) 

♦  EXP{X2}+EXP(X2}  + . +  EXP{X2)  (nA  terms) 


♦  EXP(Xk)  +  EXP{Xk}  + . +EXP(Xk}  (nk  terms) 

k 

This  is  the  convolution  of  £  n.  exponential  random 

1=1  1 

variables  where  there  are  n^  identical  exponential  random 
variables  having  the  failure  rate  X^ . 

The  moment  generating  function  technique  is  not  useful  in 
this  situation,  since  there  is  a  huga  amount  of  complexity 
involved.  This  section  deals  with  this  notation  using  the 
convolution  formula. 

1 •  Reliability  s hoctha  nd  Notation 

EXP{Xi}+EXP{Xi}*EXP  [X2) 

He  know  the  survival  function  for  EXP  £X 2 )  +EXP  { X2 }  and 
also  for  EXP{X  }+SXP{X).  The  next  two  subsections  present 
different  ways  to  reach  the  survival  function  for  the  short¬ 
hand  notation  above. 
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I 
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(a)  . 

Let  T1,T2,T3  be  ran dom  variables  distributed  as 
EXP(Xi),  EXP {Xi }  ,  EXP(X2}  respectively. 

Let  T 1*  *T1+T2.  This  random  variable  has  the  Erlang  (2,  A  } 
distr ibution. 

The  convolution  formula  for  T=T1*+T3  is; 

t 

FT(t)=PT/(t)  ♦  J  PTi(t-s)  fT  (s>  ds  ,  t>0. 

0  5t  1 

fT(t)«  (1  +  Xit)  e'Xlt+  |e-X^t"S^X12s  a"X]%s 

0 


F  (t)*e  {<X22“2  XxX2)  /  (X2-X  i)  2  ♦  X1X2/(X2-Xi)  t} 

-X2t 

♦«  {  Xx*/(X  i-X2  )  2  }  ,t>0. 


(b). 


Let  T1#T2,T3  be  random  variables  distributed  as 
EXP  C Xx }  ,  EXP fX  i }  #  EXP  (X  2}  respectively. 


Let  T2*=T2+T3.  This  random  variable  has  the  survival 
function 

-  ~Xlt  -X  t 

Ft  (t)  =X2/(X2-Xi)e  ♦  Xi/(Xx-X2)  eAzt  ,t20. 

'2 

The  convolution  formula  for  T=T1  +  T2'  is 
_  _  t  _ 

FT(t)  »  (t)  ♦  |  (t-s)  fy2 , ( s)  ds  ,  b>0 . 

0 
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The  result  from  the  above  is 


PT(t)»  (aA*ate)€^ltM4*%lt)  e  ,  t>0. 
where  a*=  (X^-OX^X! )  /(X2~  Xx)  3  ,  XxX2  Z/(X2  “  h  >  2  • 

®2 (A.l^“3Xj^X2)/(Xl  “  X  2  )  ^  »  ^2*2“ Xi  “^X2/(Xx”X2)^  • 

It  is  important  to  note  that  the  number  of  exponential 
terms  equals  the  number  of  dissimilar  failure  rates  and  each 
exponential  term  has  a  polynomial  coefficient  with  the 
degree  of  the  polynomial  equal  to  the  number  of  identical 
random  variables  having  the  corresponding  failure  rate. 

The  next  section  deals  with  the  convolution  of  exponen¬ 
tially  distributed  random  variables  ising  the  fact  illus¬ 
trated  above. 

3.  Introduction  of  an  algorithm  for  the  Convolution  of 
Exponen^IaI~Lives~ 

From  Subsection  2,  we  can  infer  the  fora  of  the  sur¬ 
vival  function  which  we  sought  at  the  beginning  of 
Section  C. 

SHORTHAND  NOTATION  : 


EXP  {X  ]}  +EXP  {X;  }  +  ... 

...  +EXP(X!  ) 

<ni 

terms) 

+  EXP{X2}  +  EXP{X2>  +  ... 

...  +EXP (X2» 

(Hi 

terms) 

♦  EXP{Xk}+EXP{Xk}*... 

•  • 

...  +5XP{Xk} 

<nw 

tens) 
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SURVIVAL  FUNCTION 


P(t)  =  A !  (t)  eXzt  ♦  A2(t)  SXzt  A  (t)  iV  ,  t  >0 


whare  Ai  (t)  *an*ai2  t+a  t*  + 

1  3 

A  2(t)  aa2i+a  22t+a  2  3t*+ 


n  -i 

+  a2n  t  , 

V* 

♦aan  t  , 

2 


Ak(t)=akl>ak2t+at3t2  + 


+a 


kn, 


tnk‘ 


a.  Example: 

SHORTHAND  NOTATION  : 

EXP{Xi}  +  EXP  (X2  } 

♦  EXP(X2} 

♦  EXP  {X3  }  +EXP  {X3  }♦  EXP  {X3  } 

SURVIVAL  FUNCTION  : 

P(t)*(au+ai2t)eXlt  +a21  eXzt  +  (a3i +a32  t+a33 1*)  eXzt,t>0. 
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b.  An  Algorithm  to  Compute  the  Coefficients 

The  algorithm  represented  below  develops  the 
survival  function  by  adding  one  random  variable  in  each  run 
As  an  example  ,  in  order  to  compute  the  survival  function 
for  the  convolution  of  ten  exponentially  distributed  random 
variables,  the  algorithm  is  supposed  to  run  ten  times. 

The  notation  used  in  the  algorithm  is: 

K  number  of  dissimilar  failure  rates 

X^  ith  type  failure  rate 

X.  failure  rate  for  the  currently  entering  life 
'e 

a.  kth  coefficient  on  the  jth  polynomial 
J  k 

n.  current  number  of  identical  lives  having 
the  ith  failure  rate 

nn^  number  of  random  variables  having  ith 
failure  rate. 

Ini  to.  al:  a  j  ^  ®  0  #  V  ^  ^  v  h  e  j=1,2,..K  ,  k=1,2,...nj 

n i *0  ,  Vi  vhBce  i= 1 , 2 , . .  K 
Input  :  X^  ,  \/ i  where  i=1,2,...K 

nn^ ,  V  i  where  is1,2,...K 
The  first  run  is  :  nn=l,  a,.,=1 


Algorithm  : 

n.  ■^♦1,  until  n.i  »nnj  . 

1e  e  e  e 

1.  Update  the  coefficients  :  a,  .for  k*2,3,..vn. 

7eK  \ 

a<t*X,a^  /(m-1)  ,  where  m«n,  -j  for  j»0, 1,  .,n. -1 
e  e  em-1  ’e  \ 

2.  Update  the  coefficient  :  ^ 

ni  6 

ai  *a1  (ailXl/(Xa-X1)*  £  ( j-1)  !  ^  a  1  •/  ( X2-X  } 

e  e  1f1g  e  j=2  e  J  e 

n^>0  n^>l 

3.  Update  the  other  coefficients:  a^  V i,k 

where  i#i  ,n.#0  for  i»1,2,.,K,  k=1,2,.,n, 
el  i 

1,ni*a,niV‘VXi)  ,vi  "h9re  1,1  ■ 

a1n  *  (xi  ain  /  *  vi  where  -*-e  and  ^>1 

e  e 

for  n*n^ -j,  j*1,2,.., n.-l 
c.  Example: 

Reliability  Shorthand  Notation  : 

EXPCXj  }  ♦  EXP  (x1  }  +EXP  (X2  }  ♦'EXP  (\2  } 

Let  T1,T2,T3,T4  be  random  variables  distributed  as 
EXP  {X 1 1  #  EXPC^},  EXP  {X2}  t  EXPCXz  }  respectively. 
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We  would  like  to  derive  a  formula  for  the  survival  func¬ 
tion  of  the  random  variable  T=T1+T2+r3+T4. 

The  use  of  algorithm: 

1 st  SON; 

nx  =  1,  ai:  =1 

At  the  end  of  this  run  we  have  only  one  random  variable, 
which  is  distributed  EXP(Aa}. 

2nd  HON; 


i  =1,  n  =2 
e  1 

ai  2 =  AiaA  A 

then 

ai2  =  Ai 

an  *a  u+0 

then 

a  U*1 

At  the  end  of 

the  2. 

run  we  have  the  survival 

function 

the  T'=T1+T2, 

where 

T1  and  T2  are  identically 

distri- 

buted  as  EXP{Aa}.  The  survival  function  is 
F<t)-(au+%2t)ex»5 
P(t)=  (It  Ait)  e’^  t>0 
Which  is  EBLANG  (2,  \i)  . 

3  rd  RON; 

i  =2,  n,  *2,  n.  *1 
e  1 
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a21  *0+aux1/  (X  j-X2  )  +  a  1!  X2/(Xi-X2)2 
*Ui/(Xi-  h)+  Xi X2/(Xi-X2)  2 
*Xi2/(Xi-X2>2 

«j2  *a12Aa/0*  -Ai)  than  a12  =  x:X2/  (Xa  -  Xi ) 

an  =(A2alI-a12)/(x2 -Xa)  then  a21  =  (X^-2XiX2) / (X2-Xx )  2 
Note  that,  here  the  coefficient  a  i2  is  the  updated  one. 

At  the  end  of  the  3.  run  we  have  the  survival  function 
for  the  random  variable  T"aTl+T2+T3,  where  T1,T2,T3  are  as 
defined  before. 

~Xi  t  -X2 1 

Sjiit)*  (aif-aijt)  e  +a2#  r  t>0 

where  a^*  (X23-2x1X2)  /  (A*  “Xx  )  2  >  a*3  X2X/<  X2  "X2  )  » 
a2r- Xi2/  (Xi  "X2  )  2  • 

Note  that  the  coefficients  are  identical  to  the  result  of 
Subsection  (b). 

«th  EON; 

ie=2,  nx-2,  n2  *2 

a 2 2  3 X 2a2i  then  a  22  =  Xi2 X2 /( Xx  “X2 )  2 

a21  *a21+ailXI/(X1  -X2)4'ai^2/(X1-X2)  2 
then  a21  -(x*  -3x22x2>/(Xi  ”X2>  3 
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aa2  -^lzXz/iXz’Xi)  theiJ  a12  =A1A22/(A2-X1  )  2 

®2  j  ”  ( X2  an”ai2^  f  ^ A2  “  A*  ^  *  ^®  ^  ~(X2  3  “  3a  XA22 )  /  (A2  “Ai  )  3 

it  the  end  of  the  4.  ran,  we  have  the  survival  function 
for  the  variable  T=T1 +T2+I3 +T4. 

PT(t)  =(a1j*-a12t)  iAlt+(a2:*a2£>  eA2t  ,t>0 

where  an  —  0.2  3-”3 AiA22)/(A2“Aj)3  »  aj2~AiA22/^A2”Ax^2  * 

a2i  =  (Ax3-3A12A2)/(A1'A2)3  »  a22  =Ai  2  A2  /  ( A2 -A  2 )  2  • 

Note  that  the  coefficients  are  identical  to  the  result  of 


Subsection  2 


III.  ESIti&filLIII  SHOR£RA  ND  APgROAIH  TO  SYSTEM  RELIABILITY 


This  section  deals  with  a  system  whose  components  have 
constant  failure  rates. 

Having  the  reliability  network  for  a  system  in  which  each 
component  has  an  exponential  life  and  knowing  the  probabili¬ 
ties  for  failures  (discussed  in  Appendix  E2)  makes  it  easy 
to  describe  the  system's  life. 

In  order  to  make  the  idea  clear,  ve  will  go  through  some 
examples. 


2  spares  (X:)  1  spare  (X3) 

SYSTEM  :  2  components  in  series  with  spares. 

DESCRIPTION  :  System  has  2  components  in  series.  Compo¬ 
nent  1  has  a  life  distributed  as  EXP{Xi  }  and  there  are  2 
identical  spares.  Component  2  has  a  life  distributed  as 
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EXP  { X2}  and  there  is  one  nan-identical  spare,  whose  life  is 
distributed  as  EXP{X3}. 

LIFE  : 


yEXP{Xi  +  X2} 

3  *i/(Xi  +  X2) 

^EXP{  Xi  +  X2}+^  /EXP  { X3  +X3  } 

1  Ai/(x1+X2)  4  X2/(Xi+X2)  /  7  Xi /(Xx  ♦Xa) 

/  \  / 

EXP  {Xi  +  X2}  t  EXP  ( Xj  +X3  }  + 

2  X2/  (Xi+  X2)  8  X3/<X!  +X3  ) 

\  .ZERO  ZERO 


2  X2/(X  !+X2) 


6  X3/(Xi  +  x3> 

kEXP^+Xi+^  •  ^EXP  [Xi  +X3  } 

5  Xi/(Xi+X3)  9  X3/(Xi  +X3) 

\  / 

nexpu1*x3  R 

10  X3/(Xx+X3) 

\  TPon 


‘ZERO 


Explanation  for  the  derivation  of  system’s  life: 

At  the  beginning  the  system  has  an  EXP  {X 1  +X2 }  life.  The 
failure  of  the  system  is  by  the  failure  of  component  1  with 
a  probability  of  Xi/(Xi+X2)  or  by  the  failure  of  component  2 
with  probability  of  X2  /(Xi  +  X2). 


2'4 


In  the  life  figure,  number  1  denotes  the  event  "failure 
of  component  1"  and  number  2  denotes  the  event  "failure  of 
component  2". 

If  event  1  c-ccurs,  component  1  is  replaced  by  one  of  the 
spares  and  system  again  functions  with  a  life  distributed  as 
since  the  exponential  distribution  has  the  memo¬ 
ryless  property  and  component  2  still  has  the  same  failure 
rate. 

If  the  event  2  occurs,  component  2  is  replaced  by  its 
spare  and  the  system  has  a  life  distributed  as  EXP{Ai  +A3 } . 

The  numbers  on  the  life  figure  correspond  to  the  transi¬ 
tions  that  can  occur.  The  probability  on  each  arc  shews  the 
conditional  probability  of  the  transition.  As  an  example 
event  3  can  occur  with  probability  of  Ai  /(Ai  +  A2)r  given  than 
event  1  has  occured  before. 

The  distribution  ZERO  defined  by  Esary[Ref.1]  and 
Repicky  [Ref. 2]  (also  defined  in  Appendix  E3)  enters  when 
life  is  exhausted. 

For  convenience  of  description,  it  is  helpful  to  define 
the  concept  of  path  used  in  this  paper.  Path  denotes  the 
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sequence  of  events  in  the  system's  life  from  the  starting 
point  to  the  point  where  the  system  is  not  functioning. 
Examples; 


Events  1  and  3  are  a  path,  which  denotes  a  sequence  of 
lives  for  the  system.  In  this  case,  the  system  has  3  expo¬ 
nentially  distributed  lives  EXP  (Xj  +  A2 )  ♦EXP  fX i  +a2  }  ♦EXP  *\2 } 
with  the  probability  of  Ax  /(  Xi  %  >  A1/(a1+A2)* 


Events  1,  4  and  8  form  another  path,  which  describes  a 

sequence  of  lives  for  the  system.  In  this  path,  the  system 
life  is  EXP{  A  x +A2  } +EXP{Ax  *a2  }+EXP{A!+ A3  }  ♦ZERO.  The  prob¬ 
ability  of  this  path  is 

[Ax/(Ax+A2  ](a2  /  (Ax  ♦A2][aj/(Xi  +A3  ij  • 


B.  THE  SOB  VITAL  FUNCTION  FOB  THE  SYSTEM 


In  this  section  we  will  deal  with  an  approach  to  obtain¬ 
ing  the  system's  survival  function. 

1 •  IX A*EL£2 

SYSTEM:  Two  components  in  parallel. 


EXP  CX!  +  X2} 


BXP  C\2} 


SXP  CAi } 


DBS CB IPTION :  Component  1  has  EXPC^l  life  and  component  2 
has  EXP(Xz}  life. 


There  are  two  failure  events  that  can  be  allowed.  The 
first  is  that  component  1  fails  at  some  time  t  and  component 
2  carries  the  system  for  the  rest  of  the  time.  The  other  one 


is  that  component  2  fails  at  some  time  during  the  mission 
and  component  1  carries  the  system  for  the  rest  of  the  time. 
There  are  more  events  such  as  no  failures  during  the  mission 
duration,  which  are  taken  care  of  by  the  representation. 


Now  we  have  two  paths 


No. of  Path  Weight 


Life 


1  Pi  2XP  {Xj +X2} +EXP  (X2  } 

2  P2=X2/(Xi+X2)  EXP{Xi+x2}  +  2XP(x1} 


Let  T  be  the  system's  time  to  failure  and  let  T1,T2,T3  be 
random  variables  exponentially  distributed  with  the  failure 
rates  X1+a2#x2»Ai  respectively. 

Than 


f  T1  ♦  T2 


T  =  { 


i T1  ♦  T3 


with  probability  pt 


with  probability  p4 


The  survival  function  can  ba  written  as 
F(t)  *  Pl  P,  (t)  ♦  p2  ?2  (t)  ,  t>0 

where  px  »  Ax  /  ( Axn2 )  ,  P2  =  Aa/ (Ai  ♦  X2  )  and  P2  (t)  ,  P2  (t)  denote 
the  survival  functions  for  the  shorthand  notations 
EXP  {Ai  +A2 }  +EXP(A2} ,  EXPCAi  ♦  X2  J'-SXPCAi  }  respectively. 
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The  survival  function  for  the  convolution  of  two  exponen¬ 
tially  distributed  random  variables  with  dissimilar  failure 
rates  is 

?it)  =  \z/ (Xz-Xi)  /(Xx-Aa  )  *Xz\>0 . 

If  we  do  the  substitutions  for  Fi  (t)  and  F2  (t)  as 
.1  =\i+\2,X2  **  A2  andAi  *  Ai+ A2  ,  A2  =  Ai  respectively, 

FjCt)  becomes 

~?l  (t)»  A2  /(-Ax)  e^Al+X^t+ (Ai +A2) /Ai  s^2t  ,  t>0, 
and  F^  (t)  becomes 

_  -( Ai+A2)t  -A2t 

(*'  *A1/(-A2)  s  +(Ai  *A2  )  /A2e  f  t>0. 

Then 

FttJ-Aj/JAi^Aa)  Fi  (t)  ♦A2/(Al+A2)  F  (t) 

F(t)=eXl  +eXl  -iXl+Xz^t  ,  t>0. 

The  result  gives  the  survival  function  that  is  well  known 
for  this  system. 


Another  example  is  the  2  out  of  3  system 


The  survival  function  is 


P(t)=p1F1  (t)+p2F2  (t)+p3?3(t)  ,  t>0. 

where 

_  -(Ai+A2+A3)t  -( A2+A3)t 

Fx  (t)  -  Qv+\£  /  (“At)  s  ♦  ^1^2  ^3)  /\  e 

—  "CAi+A2+A3 jt  -(Ai+A3)t 

P2  (")  -  C^AJ  /  (-A2)  e  ♦  (Ax^+A^  A2  « 

—  — ( Ai +A2+A  3  )t  -( Ax+A2)t 

P3  (t)  =  (Ai*-^)  /  (“A,)  e  ♦  CV5k2*Aj  A3  e 

If  we  do  the  necessary  cancellations,  we  can  get, 

f(*)  =e  (^2+^3)t+e*(^i+^3)t  +e"(Ai+A2)t  ^2e"(Ai+A2+A3)t^  t>Q 


as  desired. 

2.  General-  Procedure 


Having  the  algorithm  presented  in  Section  II,  we  car- 
treat  more  complicated  systems  under  similar  assumptions. 

The  procedure  is 

i.  Set  up  the  reliability  network  for  the  system. 

ii.  According  to  this  network,  set  up  the  system  life. 

iii.  Using  the  proper  reliability  shorthand  formulas  for 
the  related  convolutions  of  exponential  lives,  set  up  the 


survival  function 


APPENDIX  A 


A-  1 

This  section  contains  the  survival  functions  for  several 
shorthand  notations  which  were  derived  by  the  use  of  the 
approach  described  in  Section  II. 

A.  1.1  Shorthand  Notation  :  EXP{X}. 

Survival  Function:  F(t)  =  sXt  ,  t>0. 

A.  1.2.1  Shorthand  Notation:  EXP{Xi}+  EXP(X2} 

Survival  Function; 

F(t)»  X2/(A2-Xi)  £Xl  ♦X1/(Xi-X2)  eXzt  ,t>0. 

A. 1.2.2  Shorthand  Notation  :  EXP { X } +EXP  {  \  } . 

Survival  Function: 

F(t)=  (U  X  t)  e~Xt  r  t>0. 

A.  1.3.1  Shorthand  Notation  :  EXP  (Xj } +  EXP  (x2 } +EXP£x3  } 
Survival  Function: 

F(+)aau  e"Xlt*a21  eXzt  m31  eXjt  ,  t>0. 
where  an  *X2Xj/<X2  *xx )  <X 3 ” )  *  a21  *XiX3/  (Xx  -x2 )  (X,  -X2 )  , 


a3 1  =  A1A2/ (  Xi  “  A3 )  f  A 2  —  A 3 )  • 


a. 1.3.2  EXPfX  }  +EXP{  A  }  *  EXP(  A  } 

Survival  Function: 

F  (t)  =  (1+A  t  +  1/2A*t*)  eXt  ,  t>0. 

A.  1.3.3  EXP(Ai  }+EXP(A2)  +  EXP(A2  ) 

Survival  Function: 

F(t)*aueXlt+(a  +a  t)e"X2t  ,  t>0 
21  22 

whare  au=  A22/ <  X2“Ai  )  2  ,  ^1=  (Ax  2-2XiX2)  /  (Ax  -X2)  2 
a2  2®  Ai  X2/  ( X 1  “A  2  )  * 

A.  1.4.1  EXP { A  1}  +EXP (X2  }  «’EXP{A3}*EXP{A,  } 
Survival  Function: 

F(t)aaneXl  ♦a2IeX2t+a3r^3t+afreX4t  ,  t>0. 

where  a,-.  =  II  Aj/ n  (A. -A.)  vi=1»2,3,4 

11  j*M  j/M  J  1 

A.  1.4.2  EXP{ A  }+EXP(  \}  +  EXP {A  JfEXPf  A  } 
Survival  Function: 

F  (t)  *  ( 1+xt  + 1/2  x  2t2+1/6  X3t3)aXt  ,  t>0. 

A.  1.4.3  EXP(Xi}  +SXP{A2  }  ♦  EXP(A2}  +EKP  {X2  } 
Survival  Function: 
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F(t)=a11«  +(a21+a22t+a,3t  2)  e  ,  t>0. 

where  an=  a2  3/(X2-Xx  )  3  ,  ^  =1"au  * 

a22”  XxX2(Xx“2  X2)  /  (a  x“  X2)  2  »  a23“  XxX2  z/2  (a1  “a2  ^  • 

A.  1.4.4  EXP  & j  }  +EXP  (Ax }  ♦  EXP  (X2  } +EX P  [A2  } 

Survival  Function: 

—  -A2 1  -X2t 

F  (t)*  (au+al2t)  e  +  (a21«^2t)e  ,  t>0. 

where  au=  (  X23-3x22Xx  ) /(A2  “  X  x )  3,  a12*  XxX22/ (A2  -  Xx)  2  , 

a2r  (Xx3-3a?X2)/(Xx  -X2)  3»  aa^fcAiVUx-fc  )*• 

A.  1.4.5  EXP  (Xx  }+EXP{X2}«-EXP{X3}+EXP{A3  } 

Survival  Function: 

_  -Xxt  -X2t  -X  3t 

F(t)  =  a  e  *a  e  +  (a  t)  3  ,  t>0. 

U  a  31  32 

where  *  xz*3Z/(  x2- Xi )  <X3-Xx)2,  XxX»2/(Xx-X*  )  (X3  -X2  ) 
a313  XxX2/(Xi"X3)  (X2  “  X  3 )  X1X2A3L  1  /  C  X2  ~X2  )  (X  x  ” X  3 )  ‘ 

-V(xx  _x2)  (x2  ” X3 ) 3  3 


a32=XiX2X3/(x1-  x3)  {X2-X3)  • 


on  non  non  non  nooonnnnnonoooo 


This  section  introduces  a  Fortran  program  using  th 


algorithm  described  in  Section  II. 

A  PROGRAM  FOR  JHE  ALGORITHM  ONE  AT  A  TIME: 


THIS  IS  A  PROGRAM  TO  COMPOTE  THE  RELIABILITY  OF  A  SYSTEM 
WHICH  HAS  THE  RELIABILITY  SHORTHAND  NOTATION 
EXP  (LI)  ♦  .  .♦EXP  (LN) 

WHERE  THERE  IS  NO  RESTRICTION  FOR  THE  FAILURE  RATES. 
VARIABLES: 

A  (Ir J)  :  I.TH  TYPE  FAI LORE  RATE,  J.TH  COEFFICIENT 

ON  THE  POLINOH. 

NI  (I)  :  AMOONT  OF  LIVES  HAVING  THE  I.TH  TYPE 

FAILURE  RATE. 

NINIT(I)  :  AOXILAORY  ARRAY  POR  THE  NI  (I)  .  INPUT 

FOR  THE  PROGRAM. 

L  (I)  :  THE  ARRAY  FOR  THE  FAILURE  RATES 


REAL 

INTEG 


LJ20),  A  (20  *  10) 
SR  NlNIT  (2  0)  ,  N 


1(20) 


GET  INPUT 


READ  (5. 499)  T 
WRITE  (6 ,498)  T 
READ J5,  500)  K 
-  [6.1500) 


WRITE  (6 , 1 500)  K 

READ  (5.501)  (NI NIT  (I)  ,L(I)  ,1*1,  K) 
WRITS  (6 , 1 50  1)  (NINIT(I)  ,L(I)  ,1*1, K) 


2 

4 

6 

9 

8 


KK=1 


COMPUTE  COEFFICIENTS  ONE  ALL 
DISSIMILAR 


A  (KK.1)  *1. 

NI  (KK)  *1 

JJ*1 

IFJJJ.EQ.  KK)  30  TO  6 
A  iKKi  Uf *  (KK,  1)  *L  (JJ)  /(L  ( JJ)  -L  (KK)  ) 
I F  (J  J  .  EQ  .  K )  GO  TO  9 
JJ*JJ+1 
GO  TO  4 

IF  (KK.EQ.  K)  GO  TO  8 

KK-KK+1 

GO  TO  2 

CONTINUE 


BEGIN  TO  ADD  ONE  AT  A  TIME 

I  E=1 

32  IF  (NINIT(IE).  EO  .NI  (IS) )  GO  TO  99 
NI  (IE)  *NI  (IE)  +1 


UPDATE  IE 


J  =  0 

NNNN*NI  (IE)  -J 
&  JIE.NNNN)  =L  (IE)  *&  (IE.NNNN-1)  /FLOAT  (NNNN-1) 


&  (IE.NNNN)  =L  (IE)  *A  (IE. 
IF  (NNNN.  EQ.2)  GO  10  20 
J=J+1 


GO  TO  12 
20  SUM=0. 

DO  101  1=1. K 

IF(I.EQ.IE).  GO  TO  101 

S0fl*S0fl+A  (I,  II  *  L (I)  /  (Lll)  -L  (IE) ) 

IF  (NI  (I)  .  LT*  2)  GO  TO  i01 
FACT=1. 

NKK  =  NI  (I) 

DO  102  11=2, NKK 
FACT=FACT<II-1) 

SUM=SUM+A(I  ,11)  *F  ACT*L  (IS)  /  (L  (I)  -L  (IE)  )  **II 
02  CONTINUE 

01  CONTINUE 

A  (IE,  1)  =SUM+A(IE,  1) 


UPDATE  I.NE.IE 


1  =  1 

21  IF  (I. 


IF  (I.  EQ.IE)  GO  TO  25 

A  (I,  NI7I)  ) =A  (I,  NI  (I)  )  *L  (IE)  /(L  (IE)  -L  (I) 
IF  (Nlfl)  .EQ.  1)  GO  TO  26 
J=i 


NNI=NI  (I)  -J 
A  (I  .  NNl )  =  (L 


(I,NNI)  =  (L JIE^  *A  (I,  NNI)  -FLOAT  (NNI)  *A  (I,  NNI+1)  ) 
F ( J. EQ.  (n{(I)(- ij  )  GoHo  26 


IF  (  J.  EQ.  (NI  (I)  « 
J  —J  + 1 
GO  TO  24 

26  IF  (I .  EQ.  K)  GO  TO  32 
1=1+1 
GO  TO  21 
99  CONTINUE 

IF(IE.EQ.K)  GO  TO  35 

IE=IE+1 

GO  TO  32 


CALCULATE  THE  PROBABILITY 


35  CONTINUE 
PRO* 5. 

DO  40  1*1, K 
S  U  H=  0 . 

NKK*NI  (I) 

DO  4  1  J=  1 ,  NK  K 
SUM*SUM+A  (I,  J)  *T*<  J-1) 
41  CONTINUE 

PRO*PRO+SUM*EXP  (-L(I)  *T) 
40  CONTINUE 


PRINT  OUT  THE  RESULTS 


WRITE  (6,600) 

DO  17  1*1, K 
NNK=  NI  (I) 

WRITE  (6,1600)  L  (I)  ,NI  (I) 

17  SSltol 

STOP 

498  FORHAT(*  1*  ,10X, 'THE  INPUT  IS:  »  ,  //1 5X  ,  »  TIME= «  ,?  10.  4) 

499  FORHAT  F10.4) 

50:  FORMAT*  3X, 13 

501  FORMAT  3X, 13, FI  0.5) 


AMOUNT  COEFFICIENTS 


600  PORMAT ( •  1  •  ,  1 0X,  ' 

601%6^MAt'('O^ToH  ' 


FAILURE  RATE 
PROBABILITY  (  TIME  >  *,F10.5,*)=  * 


♦FI  0.7) 

6  03  FORMAT  ('O'  ,35X,  3(5  (El  2.  5,2X1  ,/)  ) 

1500  FORMAT? '  0  *  ,  14X,  '  K  = •  ,3X,I2,7X, '  (THE  NUMBER  OF 
♦DISSIMILAR  FAILURE  RATES)  f) 

1501  FORM AT ( ' 0 1 , 15X, • *  OF  R.V.  FAILURE  RATE*,/15X, 
125  (*  -*)  . //20 ( 18  X.1 3, 5X.F10.5,//) ) 

1600  FOSMATVOSlOX,-**  *  "L  -~‘ 


END 


'  A  f  i  ^  A  •  £  IV/* 

F10.5,  5X,I3) 


EXAMPLE  : 


INPUT  FOR  THE  PROGRAM 


K= 


8. 


6 

5 

5 

5 

5 

5 

5 


6.5 

1.5 

2.5 

3.5 
5.0 
9.0 


OUTPUT  FROM  5 HE  PROGRAM  : 


THE  INPUT  IS: 
TIME*  8.0000 

K  =  6 

NO.  OF  R.V. 

•  —  • 

5 

5 

5 

5 

5 

5 

FAILURE  RATE 


(THE  NUMBER  OF  DISSIMILAR  LAMDA" 
FAILURE  RATE 


6.5 
1.  5 

2.5 
3.  5 
5.  0 
9.0 
AMOUNT 


COEFFICIENTS 


S) 


6.5 

5 

0. 15713E+08 

0 .843322*07 

0. 18530E+07 

0.20038E+06 

0  .927  01E+04 

1.5 

5 

0. 27364E+  09 

-0 .89949E*03 

0. 11817E+08 

-0.73813E+06 

0 . 1 8586E+05 

2.5 

5 

-0 • 1 8295  E+ 1 1 

0 .823  12E+ 1  0 

-0.  17601E+  1  0 

0. 17341 E+ 09 

-0.1 1979E+08 

3.5 

5 

0.  18841E+  11 

0.89164E+10 

0. 25648E+ * 

0.25077E+09 

0  .  334  06E+09 

5.0 

5 

-0.83480E+09 

-0 .52324E+09 

-0.  16061 E+  09 

-0.21713E+08 

-0 .22945EI-07 

9.0 

5 

-0.39700E+05 

-0 . 15168E+05 

-0. 32323E+06 

-0.33727E+06 

-0 .27860E-01 

PROBABILITY  (  TIME  >  9.0  )=  0.6935042 


non  noonnnnnnonnnnonnnnnnnnnonnnn 


APPENDIX  B 


This  section  contains  a  program  described  in  Section  III 
to  compute  the  reliability  of  a  system. 


THIS  PROGRAM  CALCULATES  THE  RELIABILITY  FUNCTION: 


PI 

?  (T) 
* 

FI  (T)  | 

[=EXP  (LI  ,  1 ) 

+  EXP  | 

LI, 2)  + . 

.  . +EXP 

♦ 

P2 

* 

F2  T)  1 

=EXP  (L2,1 

+  EXP  I 

L2 , 2)  ♦ . . 

.  . ♦EXP 

♦ 

♦ 

P3 

• 

* 

•  • 

F3  <T)  I 

•  •  • 

[  =  EXP  (L3 , 1 ) 
•  •  •  •  • 

♦  EXP  | 

•  •  < 

[L3, 2)  ♦ . 

»••••• 

..♦EXP 

+ 

• 

• 

•  • 

•  • 

•  •  • 

•  •  • 

•  •  •  •  « 

«  •  •  •  • 

•  • 

•  •  0 

•  •  •  •  •  • 

• 

♦ 

+ 

♦ 

PI 

• 

* 

•  • 

FI  (T)  (=EXP  (LI  ,  1 )  +EXP  (LI,  2)  ♦ . 

.  . +EXP 

♦ 

• 

• 

•  • 

•  • 

•  •  • 

•  •  • 

•  •  •  •  « 

•  •  •  •  • 

m  • 

m  m  m 

• 

♦ 

PN 

* 

FN  (T)  (=EXP  (LN ,  1 )  +EXP  (LN,2)  ♦ . 

.  .  +EXP 

ill!! 


LI  ,  M 1 
L2  ‘ 
L3 


p !  I 

MI  ( 


IK  ( 
PROBA 
TIME 
A  (I*  K) 


ARRAY  FOR  ALL  LAMDAS. 
ARRAY  FOR  ALL  PROBABI 
NUMBER  OF  EXP  IN  EACH 
TOTAL  NUMBER  OF  ROWS. 
PROBABILITY  OF  SYSTEM 
TIME 

THE  COEFFICIENT  FOR  T 
I.TH  TYPE  FAILURE 
RATE  AND  K.TH  COEFFICIENT 


LITIES. 

ROW. 


SURVIVAL  AT  TIME  T. 

HE  CURRENT  PATH, 

ON  THE  THIS  POLINOMIAL 


1 


REAL  P (50) 

CALL  READ  1  (IK,P  ,T) 

S  U  M=  0 . 

DO  1  1*1, IK 
SUM=SUM+P  (I) 

IF  (ABS  (SUM-  1 . )  .  GT.  1. 3  E-5)  GOTO  199 
PROB  A=0 . 0 

DO  2  1=1, IK 
WRITE(6, 1009)  1,1.  P  (I) 

CALL  ONECO  N  (  PRO  .  T) 

2  PROBA*  PROBA+PRO*P  (  I) 

WRITE  (6,601)  T,  PRO 
STOP 

199  WRITE  (6,198)  SUM 
STOP 

601  FORMAT('1*,10X,  ’PROBABILITY  (  T 
*F1 0. 7)  ' 

10  09  FORM  AT  ('O’  ,  10  X,  12  ,  • .  '  ,  3X,  '  P  ( '  ,  I 
198  FO P.H AT ( *  1  '  ,  1 0 X,  'THE  SUM  OF  THE 
*  EQUAL  TO  1.0',/IOX,  •  SUM  *  ',F1 
END 


IMS  >  « ,F10. 5, ') =  » , 

p^o^abIlitie^  IS  NOT 
0.7) 


GET  INPUT 


SUBROUTINE  RE  AD  1  (I  K,  P  ,T) 
REAL  P  (50) 

READ  (5,  499)  T 


nnn  non  n  non  non 


WRITE(6 .1499)  T 

READ  (5,  SO  0)  IK 

WRITS  (6 ,1 500)  IK 

READ  (5,  5)  (P  (I)  ,1*1,  IK) 

RETURN 

FORM  AT (7F10 .7) 


499  FORMAT  5X,F10.4) 

500  FORMAT  5X.I5) 

1499  FORMAT  M  »,10X,  ' 

1500  FORMAT  *  0«  # 10X# » 
END 


»1  1 ,10X,  *TIME=  '  .F10.4) 
'  O'  ,  10X,  *IK  =  *3X,I3) 


SUBROUTINE  ONECON  (PRO  ,T) 
REAL  A ( 20 , 1 0).L  (2(3) 

INTEGER  NI  (20),  NINIT (  20) 
CALL  READ  (K,L, NINIT) 

CALL  ON  ED  IS  (Kf  A,NI,L) 

CALL  ONEATA(K,A,NI,L,  NINIT) 
CALL  CALPRO (K, A,NI,L,  PRO,  T) 
RETURN 
END 


SUBROUTINE  READ  (K,  L,  N  INI!) 
HEAL  L  (  20) 

INTEGER  NINIT  (20) 


:nit  (20) 


READ  (5.1)  K 
WRITE  (6.1501)  K 
READ  (5,  52)  (L  (I )  ,  NINI  T (I)  ,  1=1  , X) 

RETU  R  N 

1  F0RHAT(5X,I5) 

52  FORMAT (5X.F10.5fI5) 

1501  FORMAT  ( IQx, *  K  =*,I2,5X,*(*  OP  DISSIMILAR  FAILURE 
*  RATES.  )  •) 

END 


NINI  T (I)  ,1=1  ,K) 


SUBROUTINE  ONED  IS(  K,  A  ,  NI ,  L) 
REAL  A(20 ,10)  .L  (20) 

INTEGER  NI(20) 

KK=1 

2  AjfKK,  1)  =1. 

NI  (KK)  =  1 
JJ=1 

4  IF(JJ.EQ.KK)  GO  TO  6 

jj=jj+i 

GO  TO  4 

9  IF  (KK.EQ.K  )  GO  TO  8 
KK=KK+1 
GO  TO  2 
8  RETURN 
END 


BEGIN  TO  ADD  ONE  AT  A  TIME 


SUBROUTINE  ONEATAfK,  A  ,NI,  L,  NINIT) 
REAL  L(20)  ,A  (20.1(3) 

INTEGER  NINIT  (2  0)  ,NI(20) 

32  IF  (NINIT(IE)  .  EO.NI  (IE))  GO  TO  99 
NI(IE)=NI(IE) +T 


UPDATE 


fa 


M 


M 


12 


20 


J  =  0 

SN=N I  (12)  -  J 

«  .  vvt  if  t  L*  f 


A  (IE.NN)  *L  (IE)*AfIE,NN-1) 
IF  (N&. EQ. 2)  GO  TO  20 
J=J+1 


-1)  /FLOAT  (NN-1) 


C 

C 

C 


102 

101 


21 

24 

26 

99 

35 


GO  TO  12 
CONTINUE 

DO  101  1=  1,K 

IF  (I.  EQ  .IE)  GO  TO  101 

A(IE,1V*A  (IE » 1 )  +A(I,1)  *L(I)/(L(I) -L(IE)) 

Ir  (NI  (i)  .  LT.  2)  GO  TO  101 
FACT=T. 

NKK  =  NI(I) 

DO  102  II =2,NKK 
?  AC  T=F  AC  T  f  1 1-  1 ) 

A  (IE,  1)  *A  (IE,  1)  +  A(I  ,11)  *FACT*L  (IE)  /(L  (I)  -L  (IE)  )  **II 
1  CONTINUE 

CONTINUE 

UPDATE  I.NE.IE 

IF  (I.  20. IE)  GO  TO  26 

A  (i.  NI  (I)  )  =A(I,  NI  (I)  )  *L  (IE)  /(L  (IE)  -L  (I)  ) 
tF  (NI  (I)  .  EQ.  1)  GO  TO  26 

NNI*NI  (I)  -J 


A(X,VVX)»(£jp2)  -FLOAT  (NNI)  *A  (I,NNI+1)) 


C 

c 

c 


41 


40 


'NI  (I ) -C 

=  (i7ie)  _ 

IF  (J.EQ.  {n^i)  -  I)')”' &0  TO  26 
J*J+1 
GO  TO  24 

IF(I.EQ.K)  GO  TO  32 
1=1+  1 
GO  TO  21 
CONTINUE 

IF(IE.EQ.K)  GO  TO  35 

IE=IE+1 

GO  TO  32 

RETURN 

END 


SUBROUTINE  CALPRO  ( K,  A  ,  NI,  I,  PRO,  T) 
REAL  L(20[,  A]20 ,10) 

INTEGER  NI  (20) 

PRO=0 . 

DO  40  1*1 , K 
SUM=0. 

NKK*NI  (I) 

DO  41  J*1, NKK 
S  UM*SUN+ A  (I,  J)  *T*<  J-1) 
CONTINUE 
TTT*-L (I)  *T 
PRO=  PRO+S  0M*EXP  (TIT) 

CONTINUE 

WRITE  (6,1010)  „  ,  „ 

DO  1000  1*1,  K 


1000 

1001 

1010 


RETURN 


WRITEf^lOOl)  L ( I)  , NI  (I)  ,  (A  (I ,  J)  ,  J*1,NNN) 
CONTINUE 


FOPMAT(10X,F10.5.5X,I3,5X,  .  2  j  5  (  El  2.  5 , 5X)  ,/30X)  ) 
FORMAT ( 15xJ ’L AH dH  1 , 5  X, ' n£ * ,3 X, ' COEFr ICIENTS  ,/1 4X, 


135  ('  - ')  ) 
END 
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1 


INPOT  FOR  £HE  PROGRAM 


TIME  15 

IRONS  5 

.2 

1.ST  3 

.  1 
.8 
.4 

2.  ND  3 

.  2 
.5 
.6 

3.  RD  3 

.  1 
.5 
.3 

4.  TH  2 

.  1 
.  2 

5.  TH  3 

.  1 
.2 
.4 


2  .2 

3 
2 

1 

2 

4 
1 

3 

4 
2 

3 

3 

2 

4 
3 


2  .2 


OOTPOT  FROM  I HE  PROGRAM 


TIME*  15.0000 
IK  *  5 

1.  P  (1)  =  0.2 

K  *  3  (#  OF  DISSIMILAR 

LAM DA  HI  COEFFICIENTS 


FAILORE  RATES.) 


0.1  0  3  0.  11294E+01  0.66  342E-01  0.87075E-02 

0.80  2  0.  18742E-0  1  0.23324E-02 

0.40  1  -0.  14815E+00 

2.  P  (2)  *  0.2 

K  =3  (I  OF  DISSIMILAR  FAILORE  RATES.) 

LAM  DA  N I  COEFFICIENTS 


0.20 

0.50 

0.60 

3. 

K 

LAMDA 

oTT  o~ 

0.50 

0.30 

4. 

K 

LAMDA 


0.23  148E  +  01 
0.21 333E+02  -0, 


2  -0.  25  077E+0  2 

4  -0.  13017E+03 

1  0.  15625E  +  0 3 

P  (3)  =  0.2 

=  3  (#  OF  DISSIMILAR 

NI  COEFFICIENTS 


22222E-01 


FAILORE  RATES.) 


3  0.  7553 1E+0  1  -0.  5  4932E+00  0.27466E-01 

4  -0.  89945E+0  1  -0.9  9536E+00  -0.42847E-01 

2  0.24414E+01  -0.  1  4648E+01 

P  (  4)  *  0.  2 

*2  (♦  OF  DISSIMILAR  FAILORE  RATES.) 

NI  COEFFICIENTS 


0.10  3  0. 3  20  00E+02  -0.16  000E  +  01  3.40000E-01 

0.20  3  -0.3  10  00E+02  -0.  14000E  +  01  -0.20000E-01 

5.  P(  5)  *  0.  2 

K  *  3  (#  OF  DISSIMILAR  FAILORE  RATES.) 

LAMDA  NI  COEFFICIENTS 


0.10  2  -0.  15170E+03 

0.20  4  0.  14400E+03 

0.40  3  0. 870  37E+01 

PROBABILITY  (  TIME 


0.37  926E  +  01 

0.12  800E  +  02  0.32000E+00 
0.51 851E+00  0.  88  889E-02 
>  15. 0  ) *  0.9111188 


0 . 55556E- 01 


-0.  7  3242E-0  3 


. 10667E-01 


APPENDIX  C 


INTRODUCTION 

This  section  consists  of  a  computer  program  to  compute 
system  reliabilities  as  described  in  Section  III. 

Again  for  simplicity,  we  will  go  thruogh  an  example. 

The  structure  of  the  system  is 


Components  1,  2  and  3  have  lives  exponentially  distributed 
with  failure  rates  ,X2  ,  \3  respectively.  Also  we  have  a 
spare  for  component  1  which  is  identical  to  component  1. 

Using  the  shorthand  approach,  the  system  life  would  be 
determined  as  follows 
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Aj  /  - EXP  ^2^-3} 


EXP  {Aj+\+A3}+-A2  /(X i+A2+A£ - EXPfXi+As} 


Xi/  (x^  n3) 


EXP  {X1+X2+A3}  +W  "ZERO 


X2/ VW 


EXP  (Xj+A^  ♦' 


X 3/  (V1+x2+x3) - ZERO 


Xi/  (Xj+AjJ—EXP 


'Xs/  (Xj+Ag)— ZERO 


Some  definitions  are  necessary  before  describing  the  pro¬ 


gram 

NODE:  Each  nod9  represents  an  exponentially  distrubuted 
random  variable  with  a  certain  failure  rate.  The  number  for 
any  node  can  be  chosen  arbitrarily  but  can  not  be  used  more 
than  once. 

ABC  :  Each  arc  originates  at  a  node  and  leads  to  another 
node.  Only  one  arc  can  terminate  at  a  given  node.  Arc  num¬ 
bers  can  be  chosen  arbitrarily.  If  r.  is  the  total  number  of 
nodes  in  a  system  life  then  there  will  be  n-1  arcs  in  this 
system. 

With  these  definitions,  we  can  represent  a  system  life  in 
the  following  way 


M 


BACK  POINT 35  LIST  (I£l) :  Each  node  has  a  back  pointer.  A 
back  pointer  is  an  arc  nuiber  which  shows  which  arc  connects 
the  node  to  the  tree.  If  the  pointer  is  zero,  then  the 
related  node  is  the  root  of  the  tree. 


NODE  NO 
1 
2 

3 

4 

5 

6 

7 

8 
9 


BACK  POINTER  IPS  (I) 
0 
1 
2 

3 

4 

5 

6 

7  * 


(Here  node  1  is  the  root  of  the  tree.) 


flppE  ff)DE  LIST:  As  we  mentioned,  each  node  represents  an 

exponential  lifetime,  for  simplicity  we  can  use  some  integer 


code  numbers  for  each  failure  rate 


Code  No 
0 


Related  Fa  ilure  Ra  te-L  (I) 
(distribution  ZERO) 


1 

2 

3 


Xi  +X2  +  x3 

Xi  +  X3 
Xi  +  x2 


Here  only  the  code  number  0  is  not  arbitrary  and  0  can  be 
used  for  the  ZERO  distribution.  The  others  can  be  picked  ou4, 
arbitrarily. 


ABC  CODE  LIST:  This  list  is  similar  to  the  node  code  list. 

Arc  code  numbers  represent  probabilities. 

Code  No-  P£2& ability- P& (I) 


1  Xi/(x1n2*x3) 

2  x2/(xln2tx3 ) 

3  \3/(\i*\2*X3) 

4  X1/(x1+x3) 

5  X3/(Xj  +  X3) 


US  ORIGIN  ilSl  : 


Each  arc  has  an  origin  node  and  a  termi¬ 


nal  node 


In  the  program  we  need  to  use  only  the  origin 


Arc  No 


Origin  Node  IO(I) 


1 

2 

3 

4 

5 

6 

7 

8 

LAST  POINT  NODE  LIST  (£P)  : 
the  end  of  each  path.  Th 
list. 


1 

1 

1 

2 

2 

2 

3 

3 

This  list 
are  is  no 


indicates  the  nodes  on 
necessary  order  in  the 


LAST  22ISS  DEPTH  L£ST  (IPD)  :  This  list  indicates  the  number 
of  arcs  from  last  node  to  the  root  of  the  tree. 


L.  ?.  Depth  List  IPD (I) 


I  ?•  Node  Ljist  LP  (I) 


THIS  PROGRAM  COMPUTES  THE  RELIABILITY  CF  A  SYSTEM 
WHICH  HAS  A  SURVIVAL  FUNCTION  LIKE  BELOW: 


7(T)  = 


PI  *  FI  (T)  (*EXP  (LI  f  1)  +EXP  (L1.2)  ♦ . +  EXP(L1,M1)) 

♦  P2  *  F2(T)  (=BXP(L2,  1)  *EXP(L2,2)  ♦ . ♦  EX?(L2,M2)} 

♦  P3  *  PI  <T)  (=EXP(L3/  D  «“EXP(L3,2)  ♦ . +  EXP(L3,M3|) 


♦  PI  *  FI  (T)  (=EXP  (LI,  1 )  +EXP (LI, 2)  ♦ . +2XP  (LI  r  MI)  ) 

♦  . . . 

*■  . . . . 


♦  PN  *  FN  (T)  (=EXP  (LN  ,  1 1  +EXP (LN, 2)  ♦ . +EXP  (LN  ,  MN)  ) 


IPB  (I) 

iiWP 

I  CL  (I) 

m 

Eff0 

N 

M 

NLAST 

NDIF 

MDIF 

NPATH  () 

BPATHO 

3TM 

IP 

PROBA 

PRO 

P 

LCODE  () 

NI^T() 

M#> 

T 


VARIABLES  IN  PROGRAM: 

:  BACK  POINTER  OF  THE  I.TH  NODE 
:  PROBABILITY  CODE  OP  THE  K.TH  ART 
:  ORIGIN  NODE  OF  K.TH  ARC 
:  FAILURE  RATE  CODE  OF  THE  I.TH  NODE 
:  DIFFERENT  PROBABILITY  LIST 
:  LAST  POINT  NODE  LIST 
J  LAST  POINT  NODE  DEPTH  LIST 
:  FAILURE  RATE  LIST 
:  NUMBER  OF  NODES  IN  TREE 
:  NUMBER  OF  ARCS  IN  TREE 
:  NUMBER  OF  LAST  POINT  NODES 
:  NUMBER  OF  DIFFERENT  FAILURE  RATES 
:  NUMBER  OF  DIFFERENT  PROBABILITIES 
:  NODE  LIST  IN  A  PARTICULAR  PATH 
:  ARC  LIST  IN  A  PARTICULAR  PATH 
:  THE  NUMBER  OF  THE  CURRENT  LAST  POINT  NODE 
:  THE  DEPTH  OF  THE  CURRENT  LAST  POINT  NODE 
:  THE  CURRENT  VALUE  OF  THE  SURVIVAL  FUNCTION 
:  PROB.  OF  THE  CURRENT  PATH  SURVIVAL  AT  T 
:  THE  CURRENT  PATH  PROBABILITY 
*.  FAILURE  RATE  CODES  IN  CURRENT  PATH 
:  FAILURE  RATE  LIST  IN  CURRENT  PATH 
:  *  OF  DIS.  FAILURE  RATES  IN  CURRENT  PATH 
:  NUMBER  OF  LIFES  RELATED  TO  THE  FAILURE  RATE 
:  AUXILAURY  ARRAY  FOR  NI() 

:  COEFFICIENTS  FOR  CURRENT  PATH 
:  TIME 


MAIN  PROGRAM 


REAL  L(20)  ,PA  (2  0)  ,  LL  ( 1 1)  r  PRO,  P.  PROB  A,  T 

INTEGER  IPB  (100)  .ICL  ( 100)  ,10(1)0)  ,LP(50)  ,IPD(5  0)  , 

NPATH(II),  MPATH  (10)  ,  N  I  (11 )  ,  I,  II ,  IP  A  ( 100) 

PROB A=0 . 0 

CALL  READCLP.L, PA, IPB ,ICL, 10,1  PD, NLAST ,T, I  PA) 
WRITE  (6,1013) 

DO  1  1*1, NLAST 


I  M“LP  (I) 

IP*IPD  (I) 

NPATH  (1)*IM 

DO  9  J-1.IP 

NPATH  (J  +  l)*IO  (IPB(IM)  ) 

MPATH  (J)« IPB  (IN) 

IM*N?ATH(J+1) 

9  CONTINUE 


I  P=I  F 
NPATH  (1) 
DO  9  J*1 


nnnoo  nno  nnnn  nnnn  nnnn  nnnn  nn  nnn 


CALL  PASS (PA ,  L,  ICL , NPATH, MP ATH , IP, LL, NI , IK, ?, 

WRTT3 (6 , 19  90)  £,I,P 

CALL  ONeCON ( PRO, T,  LL,NI,IK) 

PROB  A=  PROB  A+  P*PR0 


IPA) 


1  CONTINUE 

WRITE  (6,601)  T/PROBA 
STOP 

601  FORMAT  (*  0 *,  10X,  ’PROBABILITY  (  TIME  >  ’,?10.5,»)  =  », 
♦F10.7) 

1313  FOFMAT(10X, ’BEGIN  T3  CALCULATION’ . /1 1X.20 (» -» )  ) 

19  90  FORMAT  ('0*  ,  10 X,  12,  • .  '  , 5X,  '  P  (•  ,  1 2 ,  ’  )  =  •  ,F10.  7) 

END 

GET  INPUT 

SUBROUTINE  READ  (LP,L,  PA,IPB, ICL, 10, IPD, NLAST,T ,I?A) 
REAL  L^20^PAi20) 


INTEGS 
*IPA(100) 


(100)  ,  ICL  ( 100)  ,10  (130)  ,L?  (50)  ,IPD  (50)  , 
READ  TIME  AND  #  OF  NODE 

READ  (5,  1)  T, N 


READ  (5,3) 
M=N-  1 


READ  BACK  POINTERLIST. FAILURE 
RATE  CODE  LIST 

(IPB(I)  ,  ICL (I)  ,1=1  ,  N) 


READ  ARC  ORIGIN  LIST  AND  ARC 
PROBABILITY  CODE  LIST 


READ  (5,4)  (IO  (K)  ,  I  PA  (K)  ,K=1,M) 

READ  #  OF  PATH  AND  #  OF  DIFFERENT 
FAILURE  RATE  AND  #  OF  DIFFERENT 
ARC  PROBABILITIES 

HEAD  (5,5)  NLASI,MDIF,HDIF 

READ  LAST  POINT  LIST  (LP)  AND 
LAST  POINT  DEPTH  LIST  (IPD) 

READ  (5,6)  (LP(I)  ,  IPD  (I)  ,I=1,NLAST) 

READ  FAILURE  RATES 


RETURN 
1  FORMAT 

3  FORMAT 

4  FORMAT 

5  FORMAT 

6  FORMAT 

7  FORMAT 
END 


READ  (5,7)  (L  (I)  ,1=  1.NDIF) 
READ  (5 , 7)  (PA  (I)  ,I=t,MDlF) 


5X,F10.3  ,15) 

5X^1015) 

5X;i0l5) 

5X, 3151 
5X,10I5) 
5X,5F10.  5) 


THIS  SUBROUTINE  COMPUTES  NECESSARY  LISTS  TO 
CALCULATE  THE  SURVIVAL  FUNCTION  OF  A  CURRENT  PATH 


SUBROUTINE  PASS  (PA,L, ICL, NPATH, MPATH,IP,LL, NI, IK, P, 

*1 P  A) 

REAL  PA  (201  ,LL(11)  ,L  (20) 

INTEGER  NPATH  (1 1)  ,  MP A TH  (1  0)  ,  LCODE  ( 10)  , NI  ( 1  0)  , ICL  ( 100)  , 

*P=1.(100) 

DO  1  1=1, IP 
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nnn  non  onnn  non 


1  P=P*PA (IPA  (MPATH  (I)  I  > 

I PP= IP+ 1 
11=1 

IF  (ICL  (NPATHJ1)  )  .  SQ.O)  11  =  2 
DO  2  I=II,IPP 

2  LCODE(I)  =ICL  (NPAIH  (I)) 

I  K=0 

DO  3  I=II,I?P 
1 F  (LCODE  ( I)  .  EQ.  0)  SO  TO  3 
IK=IK*1 

LL  (IK)  =L  (LCODE  (I)  ) 

LLL=  1 

jj=i+i 

IP ( JJ . GT.  IPP)  SO  TO  13 
DO  4  J=JJ,IPP 

IF(LCODE(I)  .  NE. LCODE  (J)  )  GO  TO  4 
LLL=LLL+1 
LCODE  (J)  =0 
4  CONTINUE 

13  NI  (IK)  =LLL 

LCODE  (I)  =0 

3  CONTINUE 
RETURN 

END  C 

THIS  SUBROUTINE  CONTROLS  THE 
COMPUTING  OP  THE  SUR7I7AL  FUNCTION  OP  THE  CURRENT  PATH 


SUBROUTINE  ONECON(PRD ,T,L, NINIT ,K) 

REAL  A  (1 1  ,1  r‘  - 
INTEGER  NI( 

CALL  ONEDIS 
CALL  ONEATA 
CALL  CALPRO 
RETURN 
END 

THIS  ROUTINE  COMPUTES  THE  COEFFICIENTS 
AS  EACH  IS  DISSIMILAR 


Ilf  .NINIT  ( 11) 

'K,  A,NI,L) 

K,  A ,NI ,L,  NINIT) 
K,A,NI,L, PRO, I) 


SUBROUTINE  ONEDIS  (K,  A  ,NI,  L) 
REAL  A  (  11  ,10)  .L  (11) 

INTEGER  NI  (20) 

KK=1 

2  AJKK.  1)  =1. 

NI  (KK)  =  1 
JJ*1 

4  IF  JJ  J.  EQ.  KK)  GO  TO  6 

jj=jj+i 

GO  TO  4 

9  IF  (KK.  EQ.  K  )  GO  TO  3 
KK=KK+1 
GO  TO  2 
8  RETURN 
END 


BEGIN  TO  ADD  ONE  AT  A  TIME 


32 


SUBROUTINE  ONEATA(K,  A  ,  NI,  L,  NINIT) 
REAL  LJ11J . 

integer  n: 

IE*1 

IF  (NINITJIEl.E^.NI  (IE))  GO  TO  99 


llfcii  (i  1),li(  ii) 


NI  (IE)  *NI  (IE) 


UPDATE  IS 
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onn  non 


!*■ 


NN*NI(IE)-J 

A  (IE,NN)  *L  (IE)  *A(IE.NN- 
IFfNN.  EQ.2)  GO  TO  20 


iFfNN.  EQ.2)  GO  : 
J»J+1 
GO  TO  12 
20  CO  NT  I  SO  E 

DO  101  1*1, K 


1)  /FLOAT  (SN-1) 


V  I  U  1  X—  I  m 

IPjI.EQ.IEf  GO  TO  101 

A  (IE,  if  *A  (IE ,  1 )  +A  (I.1)*L(I)/(L(I)-L(IE)  ) 

IP(Nf  (I)  .  LT.  2)  GO  TO  101 

FACT*!. 


F  AC  T  *  1 
SKK*NI  (I) 

DO  102  11  =  2, NKK 
*  FACT=FACT<II-1) 

A(IE,  1)  *A  (IE,  1)  +  A  (I  ,  II)  *FACT*L(IE)/(L(I»  -L(IE))  **II 
CONTINUE 
CONTINUE 


UPDATE  I.SE.IE 


1  =  1 

21  IF jl. EQ .IS)  GO  TO  25 


IF(I.EQ.IS)  GO  TO  25 

A  (I,  Will)  )=AJI,  NI  ( I)  )  *L  (IE)  /(L  (IE)  -L  (I)  ) 
iP(Nljfl).EQ'l)  GO  TO  26 


(NIjl)  -EQ.  1) 

NNI-NI  (I)  -J 
A  (I  ,NNI)  *  (L 


IE)  *A  (I,  NNI)  -PLOAT  (NNI)  *A  (I,  NNI+1) ) 


IF  (J.EQ.  (NI 
J*J+1 
GO  TO  24 


L  (IE)  *A  (I,  NNI 
/iL(IE)  -Lfl)  > 
I  I)  -  1)  )  GO  I 


TO  26 


26  IF(I.EQ.K)  GO  TO  32 
1=1+1 
GO  TO  21 
99  CONTINUE 

IF  (IE.EQ.  K)  GO  TO  35 
IE=IE+1 
GO  TO  32 
35  RETURN 
END 

CALCULATES  THE  PROBABILITY  FOR  THE  CURRENT  PATH 

SUBROUTINE  CALPRO (K, A , NI,  L,PRO,T) 

REAL  L  ( 1 1 )  ,  Ajll ,10) 

INTEGER  NI  ( 1  T) 

PRO=0. 

DO  40  1  =  1, K 
S  U  H=  0  . 

NKK=NI  (I) 

DO  4 1  J=1, NKK 
SUM=SUM+ A  (I,  J)  *T*<  J-1) 

41  CONTINUE 
TTT=-L(I)  *T 
ERO=  PRO+S  UM*EXP  (TTT) 

40  CONTINUE 

WRITE  (6,1011) 

DO  1000  1*1, K 


NNN*NI  jl) 
WRITE  (6,101 
1000  CONTINUE 


0)  L  (I)  ,NI  (I),  (A(I,J)  ,J=1,NNN) 


RETURN 


1010  FOP.HAT  (*0*,10X,  F10 . 7,5X,I2 , 5X,  2  (5  (El  2.  5,2X)  f /2  5X)  ) 

10 11  FORMAT (15 X, *L AH DA* ,7X  ,'NI* ,5X, » COEFFICIENTS » ,/10X, 


135  ('  -  ')  ) 
END 
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INPPT 

DAJA  FOR  IIS 

PR OGg AM: 

T.N 

NO.  1 

l 

1 

13 

1 

1 

2 

1 

3 

3 

4 

NO.  6 

5 

1 

6 

2 

7 

3 

8 

3 

9 

NO.  11 

10 

3 

1 1 

2 

12 

2 

AR.  1 

1 

1 

1 

2 

2 

1 

2 

2 

3 

AR.  6 

3 

2 

4 

3 

5 

1 

5 

2 

6 

AR.  11 

6 

2 

7 

3 

NL,  NN 

6 

3 

3 

LPI  PD 

8 

3 

9 

3 

10 

3 

1 1 

3 

12 

LP6PD 

13 

3 

LAND  A 

2.0 

.  5 

1.5 

PA*  S 

.25 

.75 

1.0 

OPT  POT  FROH  HE  PROGRAM 


BEGIN  TO  CALCOLAT ION 


1  . 

LAS  DA 


0625000 


NI  COEFFICIENTS 


1.5 

2.0 

2. 

LAUDA 


2  -  0. 80000E+02  0.  24000E+02 

2  0. 81 000E+0  2  0.  18000E+02 

P(  2)  =  0.0468750 
NI  COEFFICIENTS 


1.5 

2.0 

3. 

LAN  DA 


-0. 51200E+03  0.  96000E+02 

0. 51 300E+0  3  0.  16200E+03 

P^  31=  0.  1406250 


ni  Coefficients 


0.5 

2.0 

4. 

LAUDA 


2  -0.63578E-06  0.  11852E+01 

3  0.  10000E+01  0. 81481E+00 

P(  4)*  0.046  8750 

NI  COEFFICIENTS 


1.5 

2.0 

5. 

LAN  DA 


2  -0. 51 200E+0  3  0.  96000E+02 

3  0.51300E+03  0. 16200E+03 

P  (  51=  0.  1406250 

ni  Coefficients 


0.5 

2.0 

6. 

LAN  DA 


2 

3 


-0.63578E-06  0. 11852E+01 

0.  10000E  +  01  0.  81481E+00 

P  (  6^=  0.5625000 


ni  Coefficients 


0. 18000E+02 


0. 22222E+00 


0. 18000E«-02 


0.22222E+00 


0.5  2  0. 59259E+  00  0.  8  8889E  +  00 

2.0  2  0.40741E+00  0.22222S  +  00 

PROBABILITY  (  TINE  >  2.0  )=  0.0173515 


non  non  nonoonnonnonnnnnnnnnnnnnnnoon 


APPENDIX  D 


This  section  gives  the  program  for  the  simulation  of  a 
system's  life  mentioned  in  Section  III.  The  program  uses  a 
crude  Monte  Carlo  simulation  procedure. 


THIS  PROGRAM  SIMULATES  A  SYSTEM  HAVING 
THE  RELIABILITY  FUNCTION: 


P(T)  = 

PI  *  PI  (T 

♦  P2  *  F2(T 

♦  P3  *  P3  (T 

♦  ..... 


=  EXP(L1,  1 
*EXP  L2  ' 
=  EXP  L3 


,1  + 
,1  ♦ 
#  1  J  ♦ 


EXP  LI, 2)  + . +  EXP  (LI  ,M1 

♦  EXP  L2,2  ♦ . +EXP  L2,M2 

EXP  L3,2  ♦ . +EXP  i  L3 ,  M3 


ii 


♦ 

♦ 

♦ 

♦ 

♦ 

♦ 


PI  '*  FI(T)  (*EXP(LI,.1)  ^EXP  (LI, 2)  ♦ . ♦EXP(LI,MI)) 


PN  *  PN(T)  (*EXP(LN,  1)  ♦EXP(LN,2)  ♦ . +EXP  (L"N,  HN)  ) 


Ml 

tin 

PROB 

TIME 


ARRAY  POR  ALL  LAMDAS. 

ARRAY  POR  ALL  PROBABILITIES. 

NUMBER  OP  EXP  IN  EACH  ROW. 

FIST  NUMBER  OP  LAMDAS  IN  THIS  ROW. 

TOTAL  NUMBER  OP  ROWS. 

PROBABILITY  OP  SYSTEM  SURVIVAL  AT  TIME  T. 
TIME 


REAL  LJ500)  ,P  (5  0)  ,  PRO  B 

INTEGER  MI (50) , MT (50) , MID 0 , MIUP , I , J , N 

IX=456378 

CALL  READ (L,N, TIME, MI ,MT, P) 

CHECK  FOR  SUM  OP  P'S  EQU  AL  TO  1.0 
TOT=0. 

DO  150  1=1, N 
150  TOT=TOT+P(I) 

IF  (ABS  (TOt-1.  ).GT.  1.E-5)  STOP 

CALL  CONTRO (L,N ,TIME, MI,MT,P,IX.PROB) 

CALL  OUTPUT (L,N, TIME,  HI, MT, P, PROB) 

STOP 

END 


SUBROUTINE  READ 


SUBROUTINE  RE  AD  (L,  N,  T  IMS,  MI,MT,  P) 
REAL  L  (500)  ,P(50) 

INTEGER  MI(50) , MT(53> 


53 


1 


& 


te 


3 


« 

r~r- 

t. 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 


GET  THE  TIME 

BEAD  (5,505)  TIME 
get  the  number  of  bows 
READ  (5,501)  N 
GET  ALL  P'S 

BEAD  (5,504)  (P(I)  ,1=1  ,N) 

GET  ALL  LAMDAS  IN  THE  OB  DEB  OF  ROB  BY  ROW 

MIDO=1 

MIUP=0 

DO  100  J= 1 ,  N 

IF  (J  .  NE  •  1)  MIDO=HIUPM 

MT  (J)  =MIDO 

BEAD  (5,  520)  BI(J) 

:1IUP=MIUP  +  MI  (J) 

READ  (5,503)  <L  ( I)  ,  I  =  MIDO,  MIUP) 

100  CONTINUE 
RETURN 

501  FORMAT  (51,15) 

503  FORMAT i 5F10.5) 

504  FORMAT  5F10. 7 

505  FORMAT* 5X,F10.3) 

520  FORMATi 5X, 15) 

END 


C 

C 

C 


111 

11 


C 

c 

c 


SUBROUTINE 

SUBROUTINE  CONT RO  (L,  S  ,TIME ,  MI ,  MT  ,  P  , IX, PRO) 
REAL  L(500),P(5Cf 
INTEGER  MI (50) , MT(50) 

SUM=0 • 

DO  1  1=1, N 
NN=1 000000*P  (I) 

CALL  SIMULA(NN,L,M  1(1)  ,MT(I)  ,X, IX, TIME) 

1  S  UM=SUM+  X 
PRO=SUH/1 000000 . 

RETURN 

END 

SUBROUTINE  FOR  SIMULATION 

SUBROUTINE  SI  MULA(NN,  L,  M,  MTT,  X,  IX, TIME) 

REAL  L(500)  ,RN]50) 

INTEGER  H  ,1, J , IX 
X=0. 

MMT=  MTT+M-1 

DO  1  1  1=1, NN 
Tp,ST=0 

CALL  LEXPN  (IX, RN,M  ,16807,0) 

J  J=0 

DO  111  J*MTT , MMT 
JJ=J J+  1 

TEST=TEST  +  RN  (JJ)  /L  ( J) 

IF  (TEST. GE. TIME)  X  =  1.0  +  X 
RETURN 
END 

SUBROUTINE 

SUBROUTINE  OUTPUT  (L,  N  , TIM E,  MI,  MT ,  P,PRO) 

REAL  LJ 500)  ,P  (50f  ^ 

INTEGER  MI (50) , MT(50) 


CONTRO 


OUTPUT 
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c 

c 

c 


PRINT  OUT  ALL  THE  SYSTEM 


WRITE  (6,600) 
MIUP»MI  (1^ 


WRITE  (6 


MI  DO* 


»2, 

MT  (II 


C 

C 

C 


MIUP*ilT(I 
W RITE  (6,6 
102  CONTINUE 

WRITE  (6  ,603)  TIME,  PRO 
RETURN 


♦  MI(I) -1 

2)  P(I)  ,  (L(J)  ,J*MIDO,MIUP) 


FORMAT  STATEMENTS 

6  CO  FORM  AT  f '  1  * , 5X, •  THE  SYSTEM  IS  :•) 

6  C 1  FORMAT ( '0 ', 5X , ' F  (T)  =  *  , F10 . 7,  '  *  *,• 
1  (F10.4,3X)  ,5  (//20X,5  f?10.  4,3X1)  ) 

602  FORHAT(*0»  ,9X.»  ♦  «,FT0.7.»  *  *,»EXE 
10 •  4, 3X)  .5(//20X,5(F10 .4,3  X)  )  ) 

603  FORMATf'O f.5X,' PROBABILITY  OI 
*T=  * . Flo .5  .  *  IS  1 ,  F 10.  7) 


*T=', F10.5, 
END 


EXP  L  *  • ,5 
EXP  L  a  ' ,5 (FI 
SYSTEM  SURVIVAL  AT 


MlUl  FOR  IHE  PfiOGJW 


TIME 

♦ROW 

.2 

1.  ST 

.  1 
.8 

2.  ND 

.2 

.5 

3.  FD 


15. 

5 


.5 


.2 

.8 

.2 


.  1 
.5 

4.  TH  6 

.  1 
.2 

5.  TH  9 

.  1 
.  2 


.5 

.5 

.2 


.2 

.4 

.5 

.3 

.5 

.1 

.4 

.4 


.2 

.1 

.5 

.1 

.3 

.1 

.1 

.4 


.2 
.  1 

.5 

.  1 

.2 

.2 


OUTPUT  FROM  Jflg  PROGRAM  : 
THE  SYSTEM  IS  : 


F  (T)  = 

0.2 

* 

EXP 

L 

S 

0.  1 
0.9 

0.8 

0.4 

0.  1 

0.  1 

♦ 

0.2 

* 

EXP 

L 

S 

0.2 

0.  5 

0.2 

0.5 

0.3 

0.6 

0.5 

♦ 

0.2 

* 

EXP 

L 

— 

0. 1 
0.5 

0.5 

0.5 

0.3 

0.5 

0.  1 

0.3 

0.  1 

♦ 

0.2 

* 

EXP 

L 

S 

0.  1 

0.  2 

0.2 

0.1 

0.1 

0.  2 

♦ 

0.2 

* 

EX? 

L 

S 

0.  1 
0.2 

0.2 

0.2 

0.4 

0.4 

0.  1 
0.4 

0.2 

PROBABILITY  OF  SYSTEM  SURVIVAL  AT  T=  15.00000  IS  0.911241 
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AgpEKDIX  E 


This  section  reviews  some  notions  which  are  found  in  the 
references  for  this  paper. 

E.  1  Redundant  systems  with  exponentially  lived  compo¬ 
nents: 

In  reliabili-y  analysis,  the  term  system  is  used  to 
describe  a  set  of  components  organized  to  perform  some  mis¬ 
sion.  A  system  is  redundant  if,  in  some  fashion,  some  of  the 
components  involved  act  as  back  up  for  other  components. 

A  rcugh  definition  might  be  that  a  system  is  not  redun- 
tant  if  the  failure  of  any  one  of  its  components  causes  the 
failure  of  the  system,  and  is  redundant  if  one  or  more  of 
its  components  can  fail  without  causing  the  system  to  fail. 
Thus  redundant  systems  have  the  property  that  they  can  suf¬ 
fer  damage  through  the  failure  of  some  of  their  components 
and  still  survive.  (ESART' S  ef .  1  ]) 

E. 2  First  failure  in  a  set  cf  exponentially  lived  compo¬ 


nents: 


We  have  n  components,  each  ir.deoendent  from  the  others, 
we  want  to  compute  the  probability  that  the  j.tk  component 
fails  first. 


P(j.th  component  fails  f  ir  st)  =P  (Tj  <T^  ,  vi  *i*j) 


-A  .s 


(  -  v 

P  >Tj  ,  Vi,H#  jir  *s)  Aje'J'ds 

\ 


Cn  f.  (s)  )  ]f  (si  as 

m  1  J 


-  -AvS  -A.s 
([  n  e  1  ]A,e  J  as 

n  “  n 


o J 


1^1  Xis 

V<  I  \  >  M  I  x.  ]  e  ds 

J  (=1  1  oJ  1*1  1 


*  XJ/(ll!  Xi  )# 

E.  3  Degeneracy  at  zero  (Zero  Distribution) ; 


Let  ZEHO  be  the  name  for  the  distribution  of  a  random 
variable  that  is  degenerate  at  zero. 


If  P(T0*Q)=1,  then  we  say  that  ro  has  the  distribution 
ZERO,  or 


T  ^  ZERO 
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